PAS J: Publ. Astron. Soc. Japan , 1-??, 
(c) 2009. Astronomical Society of Japan. 



Variable Eddington Factor in a Relativistic Plane-Parallel Flow 
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Abstract 

We examine the Eddington factor in an optically thick, relativistic flow accelerating in the vertical 
direction. When the gaseous flow is radiatively accelerated and there is a velocity gradient, there also 
exists a density gradient. The comoving observer sees radiation coming from a closed surface where the 
optical depth measured from the observer is unity. Such a surface, called a one-tau photo-oval, is elongated 
in the flow direction. In general, the radiation intensity emitted by the photo-oval is non-uniform, and 
the photo-oval surface has a relative velocity with respect to the position of the comoving observer. Both 
effects introduce some degree of anisotropy in the radiation field observed in the comoving frame. As a 
result, the radiation field observed by the comoving observer becomes anisotropic, and the Eddington factor 
must deviate from the usual value of 1/3. Thus, the relativistic Eddington factor generally depends on the 
optical depth r and the velocity gradient du/dr, u being the four velocity. In the case of a plane-parallel 
vertical flow, we obtain the shape of the photo-oval and calculate the Eddington factor in the optically 
thick regime. We found that the Eddington factor / is well approximated by f{r,^) = |exp(^^). This 
relativistic variable Eddington factor can be used in various relativistic radiatively-driven flows. 

Key words: accretion, accretion disks — astrophysical jets — gamma-ray bursts — radiative transfer 
— relativity 



1. Introduction 

In the moment formalism of radiation hydrodynam- 
ics, in order to close the moment equations truncated at 
a finite order, closure relations, such as the Eddington 
approximation or more complex expressions, are usu- 
ally adopted. In a non-relativistic static atmosphere 
(Chandrasekhar 1960; Mihalas 1970; Rybicki & Lightman 
1979; Mihalas & Mihalas 1984; Shu 1991; Peraiah 2002; 
Castor 2004), for example, many researchers have used 
the Eddington approximation and variable Eddington fac- 
tors as closure relations (Milne 1921; Eddington 1926; 
Kosirev 1934; Chandrasekhar 1934; Hummer & Rybicki 
1971; Wilson et al. 1972; Tamazawa et al. 1975; Unno & 
Kondo 1976; Masaki & Unno 1978; ). In relativistically 
moving flows (Mihalas et al. 1975; Mihalas et al. 1976a, b; 
Thorne 1981; Thorne et al. 1981; Flammang 1982, 1984; 
Mihalas & Mihalas 1984; Nobih et al. 1991, 1993; Kato et 
al. 1998, 2008; Castor 2004; Mihalas & Auer 2001; Park 
2001, 2006; Takahashi 2007), the Eddington approxima- 
tion is written in the comoving frame, and then trans- 
formed to the inertial frame, if necessary (Castor 1972; 
Hsieh & Spiegel 1976; Fukue et al. 1985; Sen & Wilson 
1993; Baschek et al. 1995, 1997; Kato et al. 1998, 2008). 

As is well known, the standard Eddington approxima- 
tion is valid when the radiation field is almost isotropic. 
This conditions is violated in the optically thin regime 
or the relativistic regime with strong velocity gradients. 
Furthermore, in a relativistically accelerating flow, even 
for a comoving observer, the comoving radiation field may 
become anisotropic. In order to improve the situation we 



are confronted with, instead of the usual Eddington ap- 
proximation, one should adopt a variable Eddington fac- 
tor, which depends on the velocity gradient djS/dr as well 
as the optical depth r (Fukue 2006; Fukue & Akizuki 
2006, 2007; Akizuki & Fukue 2008; Fukue 2008; Koizumi 
& Umemura 2008). 

In this paper we semi-analytically consider and derive a 
suitable form for the variable Eddington factor for a rela- 
tivistic, plane-parallel, scattering medium, using a simple 
but physically consistent treatment of the radiation field 
in the optically thick regime. 

In the next section we describe the shape of the sur- 
face where the optical depth measured by the comoving 
observer is unity, and in section 3 we derive its expression 
when the four- velocity gradient is large. In section 4 we 
numerically calculate the comoving radiation field within 
the photo-oval, and the relativistic variable Eddington fac- 
tor. In section 5, we briefly discuss the present result and 
compare with the results of previous work. The final sec- 
tion is devoted to concluding remarks. 

2. One-Tau Photo-Oval in the Comoving Frame 

Let us suppose an optically thick, relativistic radiative 
flow, which is accelerated in the vertical {z) direction, and 
a comoving observer, who moves upward with the flow, at 
an optical depth r (figure 1). 

Sufficiently deep inside the flow, where the optical depth 
is very large, the mean free path of photons is very short. 
In such a case, within the range of the mean free path, 
for the comoving observer the flow is seen to be almost 
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Fig. 1. Schematic picture of a relativistic radiative flow in 
the vertical direction. The flow is accelerated in the vertical 
{z) direction, and has a velocity gradient. The dashed curves 
are one-tau photo-ovals observed by a comoving observer. 

uniform; the velocity gradient and resultant density gra- 
dient are negligible. Hence, the mean free path is the 
same in all directions, and the shape of the surface where 
the optical depth measured from the comoving observer is 
unity, is almost a sphere; we call it a one-tau photo-sphere 
(a dashed circle in figure 1). As a result, the comoving 
radiation field is nearly isotropic, and a usual Eddington 
approximation in the comoving frame is valid. 

On the other hand, if the velocity gradient becomes very 
large, the behavior of the radiation field may be affected 
by changes in the hydrodynamics of the fiow even on a 
lengthscale comparable to that of the photon mean free 
path. Hence, the mean free path becomes longer in the 
downstream direction than in the upstream and other di- 
rections, and the shape of the one-tau range elongates in 
the downstream direction; we call it a one-tau photo- oval 
(a dashed oval in figure 1). As a result, the comoving ra- 
diation field becomes anisotropic^ and we should modify 
the usual Eddington approximation; in these conditions 
the Eddington factor may depend on the optical depth 
and the velocity gradient. 

In addition, when the optical depth is sufficiently small 
and/or the velocity gradient is suffiently large, the mean 
free path of photons in the downstream direction become 
less than unity. Hence, the shape of the one-tau range 
is open in the downstream direction; we call it a one-tau 
photo-vessel (a dashed hemi-circle in figure 1). 

In order to obtain an appropriate form of the relativistic 
variable Eddington factor in these regimes, we thus care- 
fully treat and examine the radiation field in the comoving 
frame. 

In a previous paper (Fukue 2008, hereafter refered to 
as paper I), we have examined the one-tau photo-oval, de- 
rived the variable Eddington factor in the non-relativistic 
regime, and proved that the Eddington factor decreases 
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Fig. 2. One-tau photo-oval around a comoving observer in 
the vertical (z) one-dimensional radiative flow. The comoving 
observer is located at z = zq, where the flow four speed is 
w = Wo and the optical depth is r = tq. In the s-direction, 
that forms an angle with the downstream direction, the 
mean free path is set to be /q- 

with the velocity gradient. In paper I, however, the 
treatment was restricted to the non-relativistic regime of 
^ < 0.1, where ^ {=v /c) is the flow velocity v normalized 
to the speed of light. 

In the present study, we thus examine the behavior of 
the relativistic variable Eddington factor in the relativistic 
regime of ~ 1, under the same approach of paper I. 
Instead of however, we use the four velocity u (= 7^), 
where 7 (= 1/^1 — (3'^) is the Lorentz factor, and seek the 
appropriate form of f{r^du/dr). 



3. Shape of the One-Tau 
Optically Thick Regime 
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In this section we compute the shape of the one-tau 
photo-oval in a vertical, optically thick one-dimensional 
radiative flow under the linear approximation. 

We deflne the one-tau photo-oval as the surface where 
the optical depth measured by the comoving observer is 
unity. The situation is schematically illustrated in flgure 
2. We assume that the comoving observer in the vertical 
flow is located at 2: = 2^0, where the flow four speed is u = 
uq and the optical depth is r = tq. In the s-direction, 
that forms an angle with the downstream direction, the 
mean free path of photons is Iq. The relation among these 
quantities is 

z — zo = scosO. (1) 

The continuity equation for the stationary, one- 
dimensional relativistic flow is 



pcu = J (= const.). 



(2) 
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where p is the proper gas density, u the four velocity, and 
J the mass-flow rate per unit area. In addition, the optical 
depth r is defined by 



dr : 



(3) 



where n is the opacity, which we assume constant in the 
present analysis. 

In this paper we use a linear approximation for the flow 
field; that is to say, around the position of the comoving 
observer the flow four speed is expanded as 



du 

u = uq^ —- 

az 



{z-zo). 



(4) 



It should be noted that in paper I we expanded the flow 
speed /3, and derived several results analytically, restrict- 
ing the treatment to the non-relativistic regime. Instead, 
in the present paper we expand the four velocity and 
examine the relativistic regime, although solutions must 
be calculated numerically. 

Under these assumptions, we consider two different 
regimes, as outlined below. 

3.1. Linear Regime 

We first consider the linear regime, where the density 
is expanded as 



Po^ -r 





dp du 
du dz 
p du 
u dz 



Po- 



{z - zo) 
{z - Zo) 





(5) 



where we made use of the continuity equation (2). 

Then the optical depth Tg along the 5-direction is easily 
calculated as 



/ tvpds 
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dp du 
du dz 



cosO 



Hence, the length Iq at which 
quadratic equation. 



(6) 

1 is determined by the 



dp du 
du dz 



cosO 



ll + ^Po^o -1 = 0. 



(7) 



As is easily seen from equation (7), if there is no velocity 
gradient, npol^ = 1, and the shape of the one-tau range is 
a sphere (one-tau photo-sphere). 

The quadratic equation (7) can be easily solved, giving 



kpqIq ■■ 
where 



1- Vl- 2a cos l9 
a cos 



1 du 
u dr 



d 

—- mu 

dr 



(8) 



(9) 



In figure 3 the solutions (8) of the quadratic equation 
(7) are shown for several values of the parameters a (0.02, 




Fig. 3. Shapes of one-tau photo-ovals in the hnear regime 
for several values of the parameters a (0.02, 0.1, 0.2, 0.4, 0.5), 
from spherical to elongated. 

0. 1, 0.2, 0.4, 0.5), from spherical to elongated. For a > 0.5, 
in the downstream direction there is no solution and the 
photo-oval is open, becoming a photo- vessel. It should be 
noted that our figure 3 is equivalent to figure 3 in paper 

1, although the definition of the parameter a is different. 
From the solution (8), the breakup condition, where the 

photo-oval is open in the downstream direction, is given 
by 1 — 2acos^ < 0, or 



1 du 
u dr 



1 



(10) 



Using the continuity equation (2) and the optical depth 
(3), the breakup condition (10) is also written as 



du 
dz 



cos 6 > 



nJ 
2c ' 



(11) 



We have shown that both the shape of the one-tau 
photo-oval and its breakup condition depend on the ve- 
locity gradient 

3.2. Semi-Linear Regime 

We next consider the semi-linear regime, where the den- 
sity is given by continuity equation (2). When the flow 
speed is high and the density decreases, the mean free 
path of photons becomes large, and the one-tau range ex- 
pands. In such a case, even if the flow field around the 
comoving observer is assumed to be linear in the manner 
of equation (4), there is no guarantee that, in this regime, 
the density distribution can be approximated by equation 
(5). 

Using continuity equation (2), the optical depth Tq 
along the s-direction is expressed as 



/■'» kJ 1 

Ts = tvpas = — / —as. 
Jo c Jq u 

From equations (1) and (4), we have 



du 
dz 



scosO, 



(12) 



(13) 
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Fig. 4. Shapes of one-tau photo-ovals in the semi-hnear 
regime for several values of the parameters a (0.1, 1, 2, 3, 
4, 5), from spherical to elongated 
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dz 



cosO ds^ 



and the integral is transformed as 



1 



du 
dz 



du 



(14) 



(15) 



where we use the assumption that the velocity gradient is 
constant around z = zq. This equation (15) is analytically 
integrated to give 



1 



du 
dz 



- log — 

cos6' ^0 



Uo 







du 
dr 



log — . 

cos6> ^0 



(16) 







Thus, the length /q at which = 1 is finally determined, 
with the help of equation (13) and the definition of the 
optical depth (3), by the following equations: 



= 



du 



dr 

U = Uq - 



COSU -\- UoLOg , 

^0 

cosO hipolo. 





dr 

We finally obtain 

1 du 
e ^dr 

udr 



cosO 



(17) 
(18) 

(19) 



In figure 4 the one-tau length tzpolo (19) is plotted 
as a function of for several values of the parameters 
a (= —dlnu/dr\o) (0.1, 1, 2, 3, 4, 5), from spherical 
to elongated. When the velocity gradient is small, the 
photo-oval is spherical. As the (logarithmic) velocity gra- 
dient increases, it elongates in the downstream direction. 
However, in contrast to the linear case in the previous 
subsection or in paper I, for any values of a equation (19) 
has always solutions, and the photo-oval never break up 
in the downstream direction. 




Fig. 5. Radiation field in the one-tau photo-oval around a 
comoving observer in the vertical one-dimensional radiative 
fiow. The radiation field may become anisotropic, since in 
general the emitted intensity is not uniform and it is red- 
shifted due to the velocity difference. 

4. Comoving Radiation Field and the Relativistic 
Variable Eddington Factor 

We now derive an expression for the Eddington factor 
in the comoving frame assuming that electron scattering 
is the dominant radiative process; the radiation seen by a 
comoving observer is that emitted by the one-tau photo- 
oval around him (figure 5). 

In a static and optically thick atmosphere, the radia- 
tion field is isotropic and uniform. In the present mov- 
ing atmosphere, on the other hand, there are two reasons 
for which the radiation field is not isotropic: the emit- 
ted intensity is not uniform, and it is redshifted because 
of velocity gradients. Neglecting the local emissivity, the 
frequency-integrated intensity /co observed in the comov- 
ing frame at the position zq of the comoving observer is 
related to the frequency-integrated intensity / emitted at 
the one-tau photo-oval and to the redshift z by 



(1 + 4 



4 • 



(20) 



First, the radiative intensity / in the comoving frame 
emitted from the one-tau photo-oval walls is not known 
without solving the full transfer problem. Clearly, it is 
not generally constant, but is a function of z (or r) and fi 
(= cos 6). If the optical depth is sufficiently large, the /i- 
dependence may be safely ignored, and the emitted inten- 
sity can be assumed to be isotropic. In general, however, 
the radiative intensity increases with the optical depth, 
and therefore, the intensity from the upstream direction 
is slightly larger than that from the downstream direction 
in the present one-dimensional vertical fiow. This effect 
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of non-uniformity of intensity generally acts as a force to 
accelerate a comoving observer. As long as the intensity 
is a slowly increasing function of the optical depth, its 
non-uniformity does not affect significantly the radiation 
energy density and the radiation pressure, and therefore, 
the Eddington factor (see paper I). Hence, in order to iso- 
late the effects induced by the relative velocity between 
the comoving observer and the one-tau photo-oval, in this 
paper we assume that the intensity of the radiation emit- 
ted from the one-tau photo-oval surface is constant; I = Iq. 

Second, the intensity observed by the comoving ob- 
server is redshifted (Doppler shifted) due to the velocity 
difference between the comoving observer and the one- 
tau photo-oval walls. In an accelerating fiow, where the 
fiow speed increases toward the downstream direction, the 
relative velocity is generally positive (figure 5) except for 
some special direction (^ = 7r/2). Hence, the Doppler shift 
of intensity also causes anisotropy of the radiation field at 
the position of the comoving observer. 

The relative speed between the comoving observer 
and the one-tau photo-oval walls is given by the relativistic 
summation law as 



A/3: 

where 
(3 = 

(3o = 



/?-/?o 



2' 



(21) 

(22) 
(23) 



In addition, using equations (18) and (19), the four ve- 
locity at z is expressed as 



u = UqC 



1 du 



cosO 



(24) 



Hence, all the necessary quantities are expressed as func- 
tions of uo, du/dr\o, and cos6. 

Using the relative speed above, the redshift z is ex- 
pressed as 



1 



(25) 



and the observed intensity /co is given by equation (20). 
Furthermore, the radiation energy density Eco^ the radia- 
tive fiux Fco, and the radiation pressure Pco measured by 
the comoving observer are calculated respectively as 



cEco = J 




(26) 




/co cosOdftco, 


(27) 


cPco = J 


I CO cos^ OdVtco' 


(28) 


Finally, the 


Eddinton factor in the comoving 


frame is 


given by 
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Fig. 6. Contour maps of the relativistic variable Eddington 
factor /(li, —du/dr) multiplied by 3, plotted as a function 
of the four velocity uq {x axis) and the four-velocity gradient 
—du/dT\Q (y axis). The three panels refer to different intervals 
of values of uq and — du/dr \o. 
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Fig. 7. Contour maps of the relativistic variable Eddington 
factor f(u,—du/dr) multiplied by 3, plotted as a function of 
the four velocity uq (x axis) and the logarithmic four- velocity 
gradient —d\nu/dr\o (y axis). In the relativistic regime, the 
Eddington factor does depend only on the logarithmic velocity 
gradient. 

The numerical results of this calculation are shown in 
figure 6. In the non-relativistic regime {u < 1) and for 
small velocity gradients, the Eddington factor is roughly 
constant, as expected. That is, when the velocity gra- 
dient is relatively small, the Eddington factor is almost 
1/3. As the velocity gradient becomes large, however, the 
Eddington factor decreases, as found in paper I. 

On the other hand, the behavior in the highly relativis- 
tic regime is somewhat different. As can be seen in figure 
6, in the relativistic regime, the Eddington factor assumes 
constant values on the line —du/dr oc u. Considering this 
behavior, in figure 7 we then plot 3/ as functions of the 
four velocity uq (abscissa) and the logarithmic velocity 
gradient a (= —d\n.u/dT\Q) (ordinate). As is clearly seen 
in figure 7, except for the non-relativistic regime (u < 1), 
the Eddington factor does depend only on the logarithmic 
velocity gradient a. We found that this dependence of the 
Eddington factor on a is well approximated by 



. 1 -a 1 f^du 
f = -e " = -exp -—- 
3 3 \u dr 



(30) 



In summary, the relativistic Eddington factor in the verti- 
cal fiow varies exponentially with the logarithmic velocity 
gradient in the relativistic regime. 

In the present analysis (and paper I) we have assumed 
that the comoving intensity is uniform. Hence, the above 
result is only due to the redshift effect. 

5. Discussion 

Moment equations for relativistic radiation transfer 
have been studied extensively in literature (Lindquist 
1966; Anderson & Spiegel 1972; Hsieh & Spiegel 1976; 
Thorne 1981; Thorne et al. 1981; Udey & Israel 1982; 
Schweizer 1982; Flammang 1982, 1984; Nobih et al. 1991, 
1993; Park 2001, 2006 for spherically symmetric problems; 
Takahashi 2007 for the Kerr metric). As the moment ex- 



pansion gives rise to an infinite set of equations, one must 
truncate it at a finite order by adopting a suitable closure 
relation, in order to make the transfer problem solvable. 
When truncating at the second order, for example, in an 
optically thick fiow the Eddington approximation has been 
usually adopted as a closure relation. 

Such radiation moment formalism is quite convenient, 
and it is a powerful tool for tackling problems of rela- 
tivistic radiation hydrodynamics. However, its validity is 
never known unless fully angle-dependent radiation trans- 
fer equation is solved. 

Indeed, in the relativistic regime, where the velocity 
gradient becomes large and Doppler and aberration ef- 
fects become important, even in the optically thick regime 
the isotropy of the radiation field may break down. In 
such a situation the Eddington approximation can not 
be adopted in the comoving frame (Fukue 2005; see also 
Fukue 2008). 

Instead of the usual Eddington approximation, we can 
adopt a variable Eddington factor^ which depends on the 
fiow velocity v and the velocity gradient dvjdr as well 
as the optical depth r (Fukue 2006; Fukue & Akizuki 
2006, 2007; Akizuki & Fukue 2008; Koizumi & Umemura 
2008). In one-dimensional fiows the variable Eddington 
factor fir^dv/dr) is generally defined as 



Pco = f{r,dv/dT)E,, 



(31) 



where Eco and Pco are the radiation energy density and 
the radiation pressure in the comoving frame, respectively. 
In recent studies (Fukue 2006; Fukue & Akizuki 2006, 
2007; Akizuki & Fukue 2008; Koizumi & Umemura 2008), 
the dependence of the variable Eddington factor on the 
fiow velocity has been considered. 

Fukue (2008) has firstly examined the dependence of 
the variable Eddington factor on the velocity gradient in 
the non-relativistic regime /3 < 0.1, obtaining: 



1 / 16d/? 
3 I ^15dT 



(32) 



In this paper, using a similar approach, we have nu- 
merically obtained the Eddington factor in the relativistic 
regime: 



^ . du\ 1 / 1 du 
dr J 3 \u dr 



(33) 



Unfortunately, this relation (33) holds only in the rela- 
tivistic regime and does not reduce to the relation (32) in 
the low speed limit of ^ ^ 0. However, we can easily find 
the bridging formula, that reduces to equation (32) in the 
non-relativistic limit and approaches equation (33) in the 
highly relativistic limit, as 



f- 



1 /WpdM 
3 ^'^^ \15 u dr 



(34) 



Thus, the relativistic Eddington factor in an optically 
thick, plane-parallel fiow linearly decreases as a function 
of the velocity gradient in the non-relativistic regime, and 
then exponentially decreases as a function of the logarith- 
mic velocity gradient in the highly relativistic regime. 
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6. Concluding Remarks 

In this paper, we have derived the Eddington factor in 
a plane-parahel, relativistic and opticahy thick scatter- 
ing medium that is accelerating in the vertical direction; 
we have introduced the one-tau photo-oval observed by 
the comoving observer, and then calculated the comov- 
ing radiation field and the Eddington factor. We thus 
have shown that the Eddington factor depends on the log- 
arithmic velocity gradient in the relativistic regime (30). 
Combining this result with the analytical expression of pa- 
per I, we derived a suitable expression for the Eddington 
factor (34), valid in both the non-relativistic and fully rel- 
ativistic regimes. 

In order to treat the problem in the simplest form, we 
only considered the optically thick regime, where the one- 
tau photo-oval always exists. In the optically thin regime 
or for large velocity gradients, a photo- vessel may form, 
where the optical depth becomes less than unity in the 
downstream direction. Furthermore, in order to isolate 
the effect of redshift, we have assumed that the radia- 
tive intensity emitted from the one-tau photo-oval surface 
is constant. In the optically thin regime or in the ex- 
tremely relativistic regime, the comoving intensity would 
be a quickly varing function of the optical depth. This 
would affect the behavior of the Eddington factor (cf. 
Fukue 2006; Koizumi & Umemura 2008). We will examine 
these general cases in the future. 

Under the present plane-parallel case, the velocity in- 
crement in the downstream direction causes the density 
to decrease in the same direction; this effect is essential 
for the shape of the one-tau photo-oval and vessel. In the 
spherical flow, on the other hand, there exists a geomet- 
rical dilution effect, and the situation can be somewhat 
different, as terms appear (e.g., Peraiah 2002). The 
present approach, however, can be extended to treat also 
the spherical case and we will examine it in future work. 

The treatment of the relativistic Eddington factor pre- 
sented in this paper may turn out to be applicable in 
various aspects of relativistic astrophysics with radiation 
transfer; i.e., black- hole accretion flows with supercriti- 
cal accretion rates, relativistic jets and winds driven by 
luminous central objects, relativistic explosions including 
gamma-ray bursts, neutrino transfers in supernova explo- 
sions, and various events occured in the proto universe. 

For example, optically thick, spherically symmetric 
accretion flows onto black holes were examined (e.g., 
Tamazawa et al. 1975; Begelman 1978; Vitello 1978; 
Burger & Katz 1980; Flammang 1982; Park & Miller 1991; 
Nobili et al. 1991), while spherically symmetric, relativis- 
tic radiation hydrodynamical winds have investigated by 
several researchers (e.g.. Castor 1972; Ruggles & Bath 
1979; Mihalas 1980; Quinn & Paczyhski 1985; Paczyhski 
1986; Turolla et al. 1986; Paczyhski 1990; Nobih et al. 
1994; Akizuki & Fukue 2008). Optically thick, super- 
critical accretion flow was also extensively studied under 
the one-zone steady approximation (e.g., Abramowicz et 
al. 1988; Szuszkiewicz et al. 1996; Beloborodov 1988). 
These models have been applied to neutron-star winds 
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and gamma-ray bursts, or black- hole accretion flows. As 
a closure relation, in some of these studies (e.g, Quinn & 
Paczyhski 1985; Paczyhski 1986, 1990), the diffusion ap- 
proximation in the comoving frame and LTE was assumed. 
However, the diffusion approximation may violate in the 
extreme relativistic limit with large velocity gradients. In 
some cases (e.g., Tamazawa et al. 1975; Turolla et al. 
1986; Park & Miller 1991; Nobih et al. 1994), the vari- 
able Eddington factor was used. However, the adopted 
Eddington factor did only depend on the optical depth, 
but did not depend on the velocity gradient. 

For the problem of of the black-hole accretion flows 
and winds in the supercritical accretion regime, radia- 
tion hydrodynamical equations were solved numerically 
(e.g., Eggum et al. 1985, 1988; Kley 1989; Okuda et al. 
1997; Kley & Lin 1999; Okuda & Fujita 2000; Okuda 
2002; Okuda et al. 2005; Ohsuga et al. 2005; Ohsuga 
2006). These models have been applied to relativistic 
astrophysical jets and winds. In these studies, however, 
simple assumptions have been made for the treatment of 
the radiation field; e.g., equilibrium between gas and ra- 
diation was assumed, the flux-limited diffusion approxi- 
mation was adopted, the flow velocity was on the order 
of 0.1c, or the optically thick to thin transition was not 
properly treated. Up to now, in relation to accretion-disk 
winds, no one solved the fully relativistic radiation hydro- 
dynamical equations. 

The treatment of the relativistic Eddington factor pre- 
sented in this paper will help in investigating these types 
of astrophysical problems that involve the solution of the 
fully relativistic radiation hydrodynamical equations. 
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